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Abstract 

The  general  theory  of  Ragone  plots  for  energy  storage  devices  (ESD)  is  discussed.  Ragone  plots  provide  the  available  energy  of  an 
ESD  for  constant  active  power  request.  The  qualitative  form  of  Ragone  plots  strongly  depends  on  the  type  of  storage  (battery,  capacitor, 
SMES,  flywheel,  etc.).  For  example,  the  energy  decreases  as  a  function  of  power  for  capacitive  ESD,  but  increases  for  inductive  ESD. 
Analytical  results  for  a  representative  set  of  ideal  ESD  (battery,  capacitor,  and  SMES)  are  compared.  Furthermore,  the  effect  of  leakage 
and  of  the  specific  loss  type  (Coulomb,  Stokes,  and  Newton  friction)  is  discussed  for  inductive  ESD.  Finally,  we  address  the  problem  of 
how  composite  ESD  should  be  treated,  and  illustrate  it  for  a  battery  with  inductance.  ©  2000  Elsevier  Science  S. A.  All  rights  reserved. 
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1.  Introduction 

Energy  storage  devices  (ESD)  are  characterized  by  the 
energy  and  the  power  being  available  for  a  load  [1,2].  A 
prominent  example  is  the  comparison  of  conventional  bat¬ 
teries  and  capacitors.  While  batteries  have  high  energy 
densities  (about  105J/kg  specific  energy)  but  only  low 
power  densities  (below  100  W /kg  specific  power),  capaci¬ 
tors  have  rather  high  power  densities  (about  106  W /kg) 
but  low  energy  densities  (about  100  J/kg).  Batteries, 
capacitors,  flywheels,  superconducting  magnetic  energy 
storage  devices  (SMES),  pressure  storage  devices,  etc.,  are 
thus  located  in  characteristic  regions  in  the  power-energy 
plane.  Typical  examples  are  shown  in  Fig.  1.  These  re¬ 
gions  are  related  to  specific  applications  by  energy  and 
power  requirements.  The  boundaries  of  the  regions  are 
determined  by  internal  losses  and/or  leakage,  etc.,  of  the 
various  ESD.  The  characteristic  time  of  an  application  is  of 
the  order  of  the  energy-to-power  ratio  of  the  ESD.  In  the 
log-log  plane  of  Fig.  1,  the  time  corresponds  to  straight 
lines.  Obviously,  batteries  are  useful  for  long  time  applica- 
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tions  (>  100  s),  while  conventional  capacitors  are  useful 
for  short  time  applications  (<  0.01  s). 

Since  the  efficiency  of  an  ESD  is  usually  dependent  on 
the  working  point,  a  single  device  belongs  to  a  whole 
curve  in  the  energy-power  plane  (see  inset  of  Fig.  1). 
These  so-called  Ragone  plots  ,  which  are  usually  presented 
in  a  log-log  plot,  are  standard  in  the  battery  community 
since  a  long  time  [1].  First,  they  provide  the  limit  in  the 
available  power  of  a  battery  or  a  capacitor.  Secondly,  they 
provide  the  optimum  region  of  working,  which  is  given  by 
the  part  of  the  curve  where  both  energy  and  power  are 
high.  The  aim  of  this  paper  is  to  present  a  unified  discus¬ 
sion  of  the  qualitative  behavior  of  Ragone  plots  for  differ¬ 
ent  ESD.  Here,  we  will  focus  on  the  specific  curves  rather 
than  on  the  specific  regions  where  these  curves  are  located. 
It  turns  out  also  that  the  specific  form  of  the  Ragone  curve 
depends  on  internal  loss  and  leakage  properties  of  the 
ESD.  A  typical  qualitative  behavior  of  a  Ragone  curve  is 
sketched  in  the  inset  of  Fig.  1.  Consider  for  example  a 
capacitor  or  a  battery.  The  internal  self  discharge  leads  to  a 
decrease  of  the  energy  that  can  be  utilized,  if  the  character¬ 
istic  time  of  the  application  exceeds  the  self  discharge 
time.  This  fact  corresponds  to  a  drop  of  the  Ragone  curve 
for  sufficiently  low  power.  On  the  other  hand,  the  effective 
series  resistance  leads  to  a  lower  time  limit  and  thus  to  a 
maximum  power.  It  is  clear  that,  irrespective  of  the  type  of 
ESD,  there  are  always  physical  limits  to  minimum  and 
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power  density  [W/kg] 

Fig.  1.  Ragone  plane:  available  energy  of  an  energy  storage  device  for 
fixed  power.  Different  types  of  energy  storage  devices  are  typically 
located  in  different  regions.  Characteristic  times  correspond  to  lines  with 
unity  slope.  Every  energy  storage  device  is  represented  by  a  curve  E(P) 
(inset).  Internal  dissipation  and  leakage  losses  lead  to  a  drop  of  the  energy 
for  sufficiently  high  and  low  power. 

maximum  speed  of  discharge  of  an  ESD.  These  limits  are 
reflected  in  the  low/high  power  behavior  of  the  Ragone 
plot. 

In  the  next  section,  we  introduce  the  general  class  of 
ESD  that  will  be  investigated,  and  we  propose  a  mathemat¬ 
ical  definition  of  the  Ragone  plot.  In  Section  3,  we  discuss 
two  specific  cases  of  potential  ESD,  which  may  be  inter¬ 
preted  as  ideal  battery  and  capacitor.  In  Section  4,  the 
Ragone  plot  of  a  purely  inductive  ESD  is  studied,  which 
may  be  interpreted  as  a  SMES  or  a  flywheel.  Furthermore, 
the  effect  of  various  types  of  friction  forces  in  kinetic  ESD 
is  addressed.  Section  5  provides  a  brief  discussion  of  the 
stability  problem  of  circuits  containing  a  constant  power 
load;  as  a  particular  example,  we  discuss  the  battery  with  a 
series  inductance. 


2.  Ragone  plots 

Consider  the  general  circuit  of  Fig.  2.  For  example,  the 
ESD  may  consist  of  a  voltage  source,  V(Q)  ,  depending  on 
the  stored  charge  Q,  an  internal  series  resistor,  R,  and  an 
internal  inductance,  L.  Note  that  this  ESD  can  describe 
many  kinds  of  electrical  power  sources.  For  example,  a 
current  source  delivering  a  current  70  can  be  described  by 
L  =  0,  and  R,  V  — >  with  V/R  ->  IQ.  The  ESD  is  con¬ 
nected  to  a  load  which  draws  constant  (active)  power 
P  >  0.  Of  course,  in  general  such  a  load  is  not  related  to  a 
constant  resistance  (except  for  the  battery),  but  requires 
control  engineering.  We  assume  first  that  the  load  has  no 
reactive  power.  (In  Section  5,  we  show  that  in  general  an 


additional  apparent  power  must  be  considered.)  The  cur¬ 
rent  I  and  voltage  U  at  the  load  are  then  related  nonlin- 
early  by  U  =  P/I.  Provided  reasonable  initial  conditions, 
*2(0)  =  <2o  and  *2(0)  =  *2(l,  are  given,  the  electrical  dynam¬ 
ics  is  governed  by  the  following  ordinary  differential 
equation  (ODE)  for  Q(t ), 

where  the  dot  indicates  differentiation  with  respect  to  time. 
This  equation  applies  not  only  to  electrical  ESD  but  covers 
many  kinds  of  physical  systems  (mechanical,  hydraulic, 
etc.).  For  example,  identifying  L,  R,  and  V(Q)  =  AW/AQ 
with  inertia  (mass),  Stokes  friction  constant,  and  a  (negat¬ 
ing)  force,  respectively,  the  left  hand  side  of  Eq.  (1)  is 
Newton’s  equation  for  the  coordinate  Q  of  a  mechanical 
particle  in  the  potential  W.  The  right-hand  side  describes  a 
force  acting  on  the  particle  and  leading  to  an  energy 
release  with  power  P.  A  similar  equivalence  holds  for 
flywheels  where  Q  denotes  the  angle  coordinate. 

Even  without  reference  to  a  specific  physical  interpreta¬ 
tion  of  Eq.  (1),  we  can  define  the  Ragone  curve  as  follows. 
At  time  t  =  0,  the  device  contains  the  stored  energy, 
E0  =  LQl/2  +  W(Q0).  For  t  >  0,  the  load  draws  a  con¬ 
stant  power  P  such  that  Q(t )  satisfies  Eq.  (1).  It  is  clear 
that,  for  finite  E0  and  P,  the  ESD  is  able  to  supply  this 
power  only  for  a  finite  time,  say  t  =  tj  P).  A  criterion  is 
given  either  by  the  time  when  the  storage  device  is  cleared 
or  when  the  ESD  is  no  longer  able  to  deliver  the  required 
amount  of  power.  Since  the  power  P  is  time  independent, 
the  available  energy  is  E(  P)  =  Plf(  P).  The  curve  E(  P ) 
versus  P  corresponds  to  the  Ragone  plot.  This  approach  is 
in  fact  independent  of  the  specific  ESD,  which  is  the 
reason  why  we  call  it  ‘general  theory’  of  Ragone  plots. 
Note  that  our  definition  is  different  from  the  definition 
discussed  by  Pell  and  Conway  [3,4].  These  authors  con¬ 
sider  the  energy  of  the  ESD  but  not  Rl.r(  P). 


Fig.  2.  General  circuit  associated  with  Ragone  plots:  the  energy  storage 
device  (ESD)  feeds  a  load  with  constant  (active)  power  consumption  P. 
The  ESD  contains  elements  for  energy  storage.  Due  to  constant  power, 
energy  supply  occurs  only  for  a  finite  time,  t,,(  P).  The  energy  available 
for  the  load,  E(P)=  PC,  defines  a  Ragone  plot. 
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Before  we  calculate  E(P )  for  a  few  idealized  cases,  we 
mention  the  trivial  case  without  internal  loss,  R  =  0  .  It  is 
clear  that  for  any  reasonable  ESD  the  full  energy  is  then 
available  and  ty  =  E0/P,  such  that  the  Ragone  curve 
Ei  P  )  =  E0  is  constant  for  all  P  >  0.  The  reader  can  easily 
convince  himself  that  the  results  below  will  reproduce  this 
behavior  in  the  limit  of  vanishing  losses  and  leakage. 


3.  Storage  of  potential  energy 


In  this  section,  we  consider  ESD  without  inductance 
(L  =  0  ).  We  focus  on  the  particular  cases  of  an  ideal 
battery  and  an  ideal  capacitor.  ‘Ideal'  means  that  there  is 
neither  frequency  dependence  nor  an  intrinsic  nonlinearity 
(e.g.  faradaic  contributions,  pseudo-capacitance,  etc.).  Bat¬ 
tery  and  capacitor  differ  in  their  charge  dependence  of 
V(Q)  in  Eq.  (1). 


3.1.  Battery 


The  ideal  battery  with  capacity  Q0  (inset  Fig.  3)  is 
characterized  here  by  a  constant  (charge  independent)  re¬ 
versible  cell  voltage,  V  =  UQ  if  Qq  >  Q  >  0  and  V  =  0  if 
<2  =  0.  In  a  first  step,  we  disregard  the  leakage  resistance 
Rl.  Eq.  (1)  reads  P=  UI  =  (770  —  RI)I,  where  U  is  the 
terminal  voltage  and  I  =  Q  is  the  current.  The  solutions  of 
this  quadratic  equation  are 


7,  — - 1- 

*  2R~ 


U  2 

uo 

4 It2 


(2) 


In  the  limit  P  —>  0,  the  two  branches  correspond  to  a 
discharge  current  7+  — >  U0/R  and  to  l_—>  0.  For  the  ideal 
battery,  the  constant  power  sink  can  also  be  parameterized 
by  constant  load  resistance,  PLoad.  The  two  limits  belong 


then  to  R, 


0  (short  circuit)  and  R, 


(open 


switch),  respectively.  Clearly,  in  the  context  of  the  Ragone 
plot,  we  are  interested  in  the  latter  limit,  such  that  we  have 
to  take  the  branch  with  the  minus  sign,  I  =  /_,  in  Eq.  (2). 

Now,  the  battery  is  empty  at  time  t.y_  =  Q0/I,  where  the 
initial  charge  Q0  is  related  to  the  initial  energy  by  E0  = 
QqUq.  It  is  now  easy  to  include  the  presence  of  an  ohmic 
leakage  current  into  the  discussion.  The  leakage  resistance 
Rl  increases  the  discharge  current  I  by  U0/RL.  The 
energy  being  available  for  the  load  becomes 

2  RQnP 

Eb(P)=Ptoo  =  - - .  °  ■  (3) 


Un 


Uq  -  4RP  +2UqR/Rl 


Eq.(3)  corresponds  to  the  Ragone  curve  of  the  ideal 
battery.  In  the  presence  of  leakage,  Eb(0)  =  0  and  there 
exists  a  maximum  at  P~  Ub  /  ]J2RRL  .  Without  leakage 
(R/Rl  —>  0),  the  maximum  energy  is  available  for  vanish¬ 
ingly  low  power,  Eh(  P  — »  0)  =  E0.  From  Eq.  (3),  one 
concludes  that  there  is  a  maximum  power,  Pmax  =  I/02 / 4  R, 
associated  with  an  energy  £0/2  (here,  we  neglected  a 


Fig.  3.  Solid  curves:  Ragone  plots  of  an  ideal  battery  with  and  without 
leakage  (resistance  RL).  Dashed  curve:  secondary  branch  of  the  energy- 
power  relation,  not  useful  for  a  Ragone  curve  (see  text).  Inset:  constant 
power  load  P  connected  to  a  battery  with  capacity  Q0,  ESR  R,  and 
leakage  resistance  RL. 


small  correction  due  to  leakage).  This  point  is  the  endpoint 
of  the  Ragone  curve  of  the  ideal  battery,  where  only  half 
of  the  energy  is  available  while  the  other  half  is  lost  at  the 
internal  resistance. 

Let  us  finally  express  the  Ragone  plot  for  the  battery  in 
the  dimensionless  units  eb  =  Eb/Q0U0  and  p  =  ARP  /Uq 

eb(p)  “  2  +2  R/Rl  ‘  ^ 

Ragone  curves  (Eq.  (4))  with  and  without  leakage  are 
shown  in  Fig.  3.  The  branch  belonging  to  1+  is  plotted  by 
the  dashed  curve.  A  more  detailed  description  of  batteries 
including  Tafel  polarization,  concentration  polarization, 
etc.,  is  discussed  in  a  different  way  by  Pell  and  Conway 
[3]. 

3.2.  Capacitor 

The  case  of  an  ideal  electric  capacitor  (inset  Fig.  4)  is  a 
little  more  laborious  than  the  case  of  the  ideal  battery, 
since  now  an  ODE  rather  than  an  algebraic  equation  has  to 
be  solved.  The  electric  potential  V(Q)  =  Q/C  depends 
linearly  on  the  charge  via  a  capacitance  C.  It  is  convenient 
to  derive  from  Eq.  (1)  an  equation  for  the  voltage  drop  at 
the  load,  U  =  P/1  =  Q/C  —  Rl  with  7  =  —Q.  Differenti¬ 
ating  both  expressions  for  U  with  respect  to  time,  replac¬ 
ing  7  in  the  first  expression  with  the  help  of  the  second 
one,  and  multiplying  by  U  leads  to  the  following  ODE  for 
U2 

I  RP\dU2  2 P 
U2)^-  ~~C' 
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Fig.  4.  Normalized  Ragone  curve  for  the  capacitor.  Inset:  constant  power 
load  P  connected  to  a  capacitor  with  capacitance  C  and  ESR  R. 


Separation  of  variables  and  integration  gives  the  solu¬ 
tion 


t(U) 


c 

/ 

RP  In 

'  U2  ' 

+  Uq--U2\ 

2  P 

l 

II2 

UQ 

) 

(6) 


It  turns  out  that  there  exists  a  turning  point  with  dt/dU 
=  0,  where  the  capacitor  is  no  longer  able  to  supply  the 
required  power  P.  One  finds  for  the  corresponding  value 
of  the  voltage  Uv  =  ]/RP  (assuming  positive  U ).  Note  that 
in  contrast  to  the  battery,  the  capacitor  is  not  empty  at 
t.y  =  l(Uj,  but  there  is  a  residual  energy  Er  =  2  RCP.  This 
relation  already  anticipates  the  existence  of  a  maximum 
power  Pmax  =  E0/2  RC.  The  time  tx  is  reached  when  the 
capacitor  is  no  longer  able  to  supply  the  required  power  P. 
In  order  to  calculate  the  Ragone  plot,  Ec  =  Pt one  has  to 
be  careful.  Indeed,  U0  in  Eq.  (6)  depends  on  P,  because 
Uc  =  U  +  RI  =  U  +  RP /U  is  the  voltage  drop  at  the  ca¬ 
pacitance.  Hence,  the  Ragone  curve  of  the  capacitor  is 


c  1 

/  RP 

,  ,  \ 

RPln  — 

U2J 

+  U2  -  RP  I , 

(7) 

U»-  — + 
0  2 

vf— ■ 

(8) 

where  the  initial  capacitor  voltage  Uc  0  is  related  to  the 
total  energy  by  EQ  =  CU(2{)/2.  Note,  as  for  the  battery,  we 
had  to  choose  Uc  0  out  of  two  branches.  However,  we  will 
not  discuss  the  secondary  branch,  which  is  irrelevant  for 
our  purposes.  As  mentioned  above,  there  exists  a  maxi¬ 
mum  power  Pmax  =  UqQ/4R  that  can  be  delivered  by  the 
capacitor  and  for  which  Ec  —>  0.  Our  results  differ  from 
those  of  Pell  and  Conway  [3],  who  find  a  finite  energy  for 
P  P,nax  •  The  reason  lies  in  the  different  definitions  of  the 
Ragone  curves,  E(P).  For  vanishing  power,  P —>  0,  the 
whole  energy  E0  is  available.  In  the  presence  of  leakage,  it 
holds  Ec(  P)  — >  0  for  P  — »  0,  as  for  the  battery.  However, 


this  case  will  not  be  discussed  here,  since  the  discharge 
time  of  a  capacitor  is  much  shorter  than  usual  leakage 
times.  In  the  dimensionless  units  ec  =  2Ec/CUq0  and 
p  =  2RCP/Eq,  the  Ragone  curve  reads 


edP)  = 


1 


—p  In 


(l  +  yjl  ~pf  ~p 

(i  +  iji -Py 


(9) 

(10) 


This  result  is  shown  in  Fig.  4. 


4.  Storage  of  kinetic  energy 

Inductive  or  kinetic  ESD  store  energy  exclusively  as 
kinetic  energy  (coming  from  the  L-term  of  Eq.  (1)),  i.e., 
V  =  0.  In  contrast  to  storage  of  potential  energy,  where 
load-free  losses  (P  —>  0)  are  related  to  a  separate  ‘leakage’ 
property,  the  kinetic  ESD  dissipates  energy  due  to  internal 
friction,  related  to  R  in  Eq.  (1).  In  practice,  inductive  ESD 
have  thus  very  small  friction  forces  or  internal  resistances. 
We  will  show  that  internal  friction  influences  kinetic  ESD 
only  at  low  power,  and  that  the  qualitative  behavior  of  the 
Ragone  plot  is  insensitive  to  the  specific  friction  among 
the  most  common  types.  At  high  power,  the  energy  turns 
out  to  be  limited  by  a  finite  bypass  resistance. 

4.1.  SMES 

In  the  following,  we  consider  a  SMES  where  energy  is 
stored  inductively,  i.e.,  in  the  magnetic  field  of  the  induc¬ 
tance  L  (inset  Fig.  5).  Nevertheless,  the  following  discus¬ 
sion  holds  also  for  other  kinds  of  inductive  storage  de- 


Fig.  5.  Normalized  Ragone  curves  for  the  inductive  ESD  with  Coulomb 
(C),  Stokes  (S).  and  Newton  (N)  friction.  The  dashed  double-dotted  curve 
corresponds  to  a  SMES  with  an  ohmic  bypass,  4R/Rb  =  0.001.  Inset: 
constant  power  load  connected  to  a  SMES  with  inductance  L,  ESR  R, 
and  bypass  resistance  Rb. 
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vices.  The  stored  energy  is  £0  =  L/02 /2,  where  L  is  the 
inductance  and  70  is  the  initial  loop  current.  At  time  t  =  0, 
the  loop  is  opened  by  inserting  a  high  resistance  Rb  and 
immediately  closed  by  connecting  it  in  series  to  the  con¬ 
stant  power  load.1  The  series  resistance  R  consists  of  all 
losses  in  contacts,  switches,  etc.,  except  the  losses  in  the 
load.  Consider  first  the  case  where  the  bypass  resistance 
can  be  disregarded,  Rb  -*  °°.  The  basic  Eq.  (1)  reduces  to  a 
first-order  ODE  for  the  current,  Ldl/dt  +  RI  =  —P/I. 
(By  the  way,  the  mechanical  interpretation  corresponds  to 
a  deceleration  of  a  moving  particle  with  velocity  7  at 
constant  power  gain.)  Multiplication  with  7  leads  to  the 
following  ODE  for  the  kinetic  energy,  Wkl n  =  LI2/ 2, 


dWkin  Wkin 

- + - +  7>  =  0, 

d  t  tt 


(11) 


load  current  7P  =  (/  +  \jl2  —  AP/Rh  )/2.  Consequently, 
the  differential  equation  for  the  SMES  current  becomes 

2  P 

LI  =  —RI - .  (14) 

/  +  y/2  -  *p/Rh 

It  is  cleai-  from  Eq.  (14)  that  there  exists  a  maximum 
power  value  Bmax  =  Rb  702/4.  Eq.  (14)  is  solved  by  first 
transformation  to  the  variable  lTkin,  and  then  separation  of 
the  variables  t  and  W.  One  ends  with  an  integral  expres¬ 
sion  for  l(  IT  )  ,  which  leads  to  if  P )  by  taking  W  = 
2 PL/ Rb.  The  result  for  the  Ragone  plot,  E  =  t/P ),  will 
be  presented  in  dimensionless  variables  in  the  following 
sub-section. 

4.2.  Dependence  on  the  type  of  friction 


where  tl  =  L/2R  is  the  energy  relaxation  time.  The 
solution  of  Eq.  (11)  with  initial  condition  lTkl  n(0)  =  E0  is 


WkJf) 


PL 
2  R  ‘ 


(12) 


In  the  inductive  case,  there  is  only  a  single  solution 
branch.  The  storage  device  is  empty  when  lTkjn(t)  =  0. 
Solving  this  for  l(  =  t/P))  yields  for  the  Ragone  curve 


PL 

E/P)=LP=—ln\ 


(13) 


The  ideal  inductive  ESD  show  a  behavior  fully  different 
from  an  ideal  device  storing  potential  energy  without 
leakage.  In  particular,  7q  — >  0  for  P  —>  0,  and  Ei  -»  E0  for 
P  — >  oo.  The  total  energy  is  available  for  fast  energy  con¬ 
sumption,  while  for  low  power  everything  is  lost  in  R. 
This  becomes  cleai-  in  the  framework  of  the  particle  pic¬ 
ture.  Drawing  energy  from  a  moving  particle  requires  the 
action  of  a  friction  force,  which  in  turn  diminishes  the 
particle’s  velocity.  If  this  friction  force  is  much  weaker 
than  internal  friction,  the  total  energy  is  dissipated  inter¬ 
nally  and  only  a  small  fraction  of  the  energy  is  available 
for  the  constant  power  load. 

In  reality  it  is  of  course  not  possible  to  get  arbitrarily 
large  power.  Cutoffs  must  exist  for  all  kinds  of  ESD  and 
for  both  limits  P  —>  0  and  P  — >  Similar  to  the  leakage 
resistance  of  a  battery,  which  becomes  important  at  P  —>  0, 
for  a  SMES  the  bypass  resistance  Rb  (e.g.,  via  the  power 
electronics)  must  be  taken  into  account  for  P  — »  'A  This 
means  that  the  load  current  7P  is  the  total  current  7  minus 
the  bypass  current  7b.  Since  the  voltage  across  the  load 
and  the  bypass  is  given  by  P/Ip  =  RbIb,  one  finds  for  the 


1  We  mention  that  the  switch  problem  is  technically  non-trivial;  how¬ 
ever,  here  we  are  mainly  interested  in  the  principal  response  of  an 
inductive  ESD. 


Three  different  types  of  friction  forces  are  often  ob¬ 
served  in  physical  systems.  First,  Coulomb  friction  is 
velocity  independent,  Ffr  =  —  sign(<2)7fc,  and  is  observed 
in  many  mechanical  systems.  Secondly,  Stokes  friction  is 
proportional  to  the  velocity,  Fft  =  —KSQ,  and  occurs  usu¬ 
ally  in  viscous  flow.  Thirdly,  Newton  friction  is  propor¬ 
tional  to  the  kinetic  energy,  Fb  =  —  K N  •  sign (Q)Q2,  and 
occurs  in  weakly  turbulent  flow.  Dissipation  in  electric 
systems  corresponds  to  Stokes  friction  with  Ks  =  R,  ac¬ 
cording  to  the  fact  that  a  current  can  be  understood  as  a 
viscous  flow  of  charge  carriers.  In  an  inductive  ESD, 
Stokes  friction  gives  the  Ragone  curve  (Eq.  (13)),  which 
can  be  expressed  in  the  dimensionless  units  es  =  2  E  J LI{2 
and  p  =  P /Rif  as 

es(p)  =pln|l  +  -)•  (15) 


In  the  presence  of  a  finite  bypass  resistance,  the  discus¬ 
sion  at  the  end  of  the  previous  subsection  leads  to 


f  1 

es(p)  =  Pj4l 


dw 


4  pR 

Kh  W  + 


2  P 


l  +  i  1 


4  Rp 

Ruw 


(16) 


The  result  is  shown  in  Fig.  5  by  the  dashed  double-dotted 
curve.  One  clearly  observes  a  sharp  drop  at  high  power 
P~RhIf. 

In  the  following,  we  compare  the  result  of  Eq.  (15)  with 
the  Ragone  plots  for  the  kinetic  systems  with  Coulomb  and 
Newton  friction. 

Replacing  RQ  by  the  Coulomb  friction  constant,  Kc, 
in  Eq.  (1)  with  V  =  0,  leads  to  the  ODE  LI  +  Kc  =  —P/I. 
Integration  of  this  ODE  and  putting  7(c)  =  0  yields  the 
Ragone  curve  EC(P )  =  Ptx.  In  dimensionless  units  ec  = 
Ec/Eq,  p  =  P/KcIq,  the  Ragone  plot  reads 

ec(p)  =2p|l -pln|l  +  -j  j-  (17) 
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An  analogous  calculation  leads  to  the  Ragone  curve 
En(P)  for  Newton  friction,  which  reads  in  the  dimension¬ 
less  units  eN  =  EN/E0  and  p  =  P/KnI^  as  follows 


«N  (P) 


2p 2/3 

3 


\^(arctan{2p  1/3/-/3~  —  1/^3"}  (18) 


1 

+  7r/  6)  +  —  In 


1  +  p 

(l+P1/3)3 


(19) 


The  three  cases  are  shown  in  Fig.  5.  Note  that  the 
power  is  scaled  with  respect  to  the  initial  internal  loss 
power,  70  Ffr(  t  =  0),  hence,  a  direct  comparison  of  efficien¬ 
cies  is  not  reasonable  from  this  scaled  figure. 


5.  General  cases 

The  previous  discussion  pretends  that  the  Ragone  curve 
of  an  arbitrary  dynamic  system  can  be  simply  derived 
analytically  or  numerically,  be  it  of  the  form  of  Eq.  (1),  be 
it  more  general.  For  example,  it  is  in  principle  straightfor¬ 
ward  to  derive  Ragone  plots  for  a  hybrid  system  consisting 
of  two  different  batteries  being  parallel,  or  battery  and 
capacitor,  etc.  In  general,  however,  things  are  far  more 
complicated.  ‘General’  means  that  both  potential  energy 
and  inductive  energy  appeal-.  As  an  illustration,  we  con¬ 
sider  the  simple  example  where  the  ideal  battery  is  coupled 
in  series  to  an  inductance  L  (inset  Fig.  6).  Using  the 
definitions  of  the  previous  sections,  Eq.  (1)  can  be  written 
as 

L/=-y(/-/+)(/-/-),  (20) 

where  the  I  +  are  given  in  Eq.  (2)  and  are  obviously  steady 
states  of  Eq.  (20).  Since  7+  >  7_  >  0,  it  turns  out  that  I+  is 


Fig.  6.  Normalized  Ragone  curves  of  the  battery  with  inductance  (inset) 
for  various  values  of  LU0 /SR2 Q0. 


linearly  stable,  while  7_  is  unstable.  Indeed,  replacement 
of  I  in  Eq.  (20)  by  the  weakly  perturbed  steady  states, 
7+  +  A/,  leads  to  an  ODE  for  the  perturbation  A 7.  Since 
A  7  is  assumed  to  be  small,  only  the  linear  terms  are 
considered.  The  resulting  ODE  is  of  the  form  A  7  =  A  +  A  7, 
where  +  indicates  the  steady  state  7+  under  considera¬ 
tion.  One  readily  finds  that  A+  is  negative  while  A_  is 
positive,  such  that  perturbations  of  7+  and  7_  are  damped 
out  and  increase  with  time,  respectively.  A  similar  behav¬ 
ior  is  valid  for  the  capacitor.  In  other  words:  in  the 
presence  of  an  inductance,  the  solution  needed  for  the 
Ragone  plot  (here:  7_)  is  not  stable  and  cannot  be  obtained 
by  using  a  load  with  U  =  P/I. 

This  forces  us  to  go  back  to  physics  and  to  the  mod¬ 
elling  of  an  element  which  draws  a  constant  active  power. 
Let  us  recall:  for  the  ideal  battery  without  inductance,  the 
constant  power  sink  is  just  a  constant  resistance  (active 
power  sink).  In  the  general  case,  however,  the  load  could 
be  a  general  impedance  with  a  certain  apparent  power.  An 
appropriate  reactance  could  remove  the  instability.  In  the 
following,  we  sketch  how  this  works  for  the  battery  with 
inductance  by  an  infinitesimal  over-compensation  of  the 
inductance. 

Consider  thus  a  load  in  Eq.  (1)  with  a  voltage 
U  =  Pq/1  —  Lp7,  instead  of  only  U  =  P /I.  We  do  not  ask 
how  the  load  is  able  to  keep  this  voltage-current  relation, 
but  consider  it  rather  as  a  black  box  containing  probably 
some  complicated  power  electronics  for  appropriate  con¬ 
trol.  Furthermore,  since  we  are  mathematically  faced  with 
a  first-order  ODE,  we  have  to  specify  the  realistic  initial 
conditions.  We  suppose  that  7(7  — »  0)  — >  0.  The  load  power 
is  given  by  P  =  UI  =  PQ  —  (Lp/2)dl2/dt.  Because  7  is  a 
function  of  time,  P  is  not  constant  as  required.  We  take 
Lp  =  L  +  8L  with  8 L  positive  and  very  small.  Hence,  L 
has  to  be  replaced  by  —  8L  in  Eq.  (20).  Due  to  the  change 
of  the  sign  of  the  coefficient  of  7,  the  steady  states  7  + 
interchange  their  stability  properties:  7_  becomes  stable. 
Due  to  the  smallness  of  8  L,  the  relaxation  time  which  is 
proportional  to  8L  can  be  made  arbitrarily  small.  Summa¬ 
rizing,  7  relaxes  very  fast  to  the  desired  solution  7_.  In  the 
limit  8 L  — >  0,  the  power  is  practically  constant,  P  =  P0, 
except  in  a  small  initial  time  interval.  Up  to  t  =  A,  the 
energy  gained  by  the  load  is  E(P)  =  Ptx  —  LI2  / 2,  were 
we  used  7(0)  =  0  and  7(a)  =  7(=  I  ).  The  time  A  is 
obtained  from  Eq.  (3).  We  recognize  that  for  constant 
current  7,  the  energy  LI2/ 2  is  stored  in  the  inductance, 
which  reduces  the  total  amount  of  available  energy.  How¬ 
ever,  at  t  =  A,  this  secondary  energy  can  be  recovered  by 
a  zero  resistance  bypass  across  the  battery,  and  running  the 
whole  device  as  an  inductive  ESD  with  initial  energy 
LI2 /2.  Using  the  result  of  Eq.  (13),  we  finally  obtain  the 
Ragone  curve  (primary  and  secondary  energy) 
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Fig.  7.  Maximum  available  battery  power  as  a  function  of  the  normalized 
inductance. 


with  I  =  I_(P )  from  Eq.  (2).  As  one  expects,  if  L  is 
small,  one  has  only  a  weak  correction  to  the  Ragone  curve 
of  the  pure  battery;  this  fact  is  not  mathematically  trivial 
since  adding  an  inductance  corresponds  to  a  singular  per¬ 
turbation  of  a  device  storing  only  potential  energy.  The 
result  of  Eq.  (21)  can  be  expressed  in  dimensionless  units, 
eb,  =  Ebi/U0QQ  and  p  =  4RP/U2. 
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which  is  shown  in  Fig.  6  for  various  values  of  the  induc¬ 
tance  in  units  of  &R2Q0/Uq.  In  general,  the  presence  of  an 
inductance  lowers  the  available  energy  and  the  maximum 
power  of  a  battery.  Moreover,  it  turns  out  that  for  suffi¬ 
ciently  large  inductance,  Ebi  =  0  for  P  <  Uq  / 4  R.  In  Fig. 
7,  the  dependence  of  Pmax  on  L  is  shown  in  dimensionless 
units. 


E2(P)...,  the  Ragone  plot  belongs  to  the  maximum  en¬ 
ergy,  E(P)  =  max  n{  En(  P)}. 

In  the  case  of  an  ESD  containing  purely  kinetic  energy 
or  purely  potential  energy,  the  power  sink  can  be  modelled 
by  a  voltage  (or  a  force)  equal  to  —P/Q,  where  Q  is  the 
velocity  of  the  dynamic  variable.  For  the  battery  and  the 
SMES,  additional  loss  mechanisms  as  leakage  and  ohmic 
bypass,  respectively,  were  included. 

It  turns  out  that,  for  realistic  cases,  a  finite  available 
energy  is  restricted  to  a  finite  power  region  0  <  P  <  Pmax . 
For  ESD  storing  potential  energy,  the  high  and  low  power 
limit  is  determined  by  internal  friction  losses  and  leakage, 
respectively.  For  ESD  storing  kinetic  energy,  the  low  and 
high  power  limit  is  determined  by  internal  friction  and 
bypass  losses,  respectively.  As  one  result,  a  consideration 
of  practical  Ragone  plots  allows  to  draw  conclusions  for 
leakage  and  loss  mechanisms.  As  another  result  important 
for  engineers,  it  is  possible  to  figure  out  bounds  for 
leakage,  loss,  etc.,  if  one  focusses  on  a  specific  applica¬ 
tion. 

In  the  general  case  of  mixed  energy  devices,  the  load  is 
a  more  complicated  device,  containing  reactive  parts,  in 
order  to  keep  a  constant  active  power.  As  an  illustrative 
example,  we  discussed  the  battery  with  series  inductance, 
where  the  load  has  to  compensate  the  inductive  voltage. 
For  the  sake  of  clearness,  the  examples  discussed  in  this 
paper  were  throughout  analytically  solvable. 

We  note  that  there  exist  different  ways  to  calculate 
energy-power  relations  of  ESD.  An  example  is  provided 
by  a  very  convenient  linear  approach  in  the  case  where  the 
information  on  the  frequency  dependent  impedance  of  the 
ESD  is  given  in  Ref.  [5].  However,  we  believe  the  present 
approach  to  be  the  most  natural  one. 
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